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PACS numbers: 74.45.-|-c,74.40.-|-k,74.72.-h,74.50.-|-r 



I. INTRODUCTION 

It is now widely accepted that the order parameter 
in the superconducting state of the cuprates has (1^2 _y2 
symmetry. Several experiments, including the tri-crystal 
ring experiments^ corner junction experiments;^ and the 
c-axis STM probes of impurity states^*^ have all pro- 
vided strong evidence for the broken reflection symmetry 
of the d^2_y2 order parameter. 

Experiments based on tunneUing within the a6-plane 
are another class of experiments which provide valuable 
information about the properties of the cuprates. Con- 
ductance measurements on planar or point contact (by 
STM) normal metal — high- Tc superconducting junctions, 
and on Josephson junctions (e.g. grain boundaries) probe 
the electronic states of the cuprates near the surface or 
at the interface. In addition to changes in the atomic- 
scale structure of the interface, the dx2-y2 superconduct- 
ing state is highly distorted by interface scattering and 
disorder. The electronic spectrum is strongly modified, 
and the dx^-y^ order parameter is in general suppressed 
on the coherence length scale. In fact, a standard fea- 
ture of the a6-plane conductance is a large zero-bias con- 
ductance peak (ZBCP)iSiL& Its sensitivity to impurity 
scattering and the splitting of the peak in a magnetic 
field agree well with theoretical predictionsSiiSiii of sur- 
face Andreev bound states with large spectral weight that 
provide a resonant channel for tunnelling near zero bias. 
Normal scattering of quasiparticles by the surface from 
a positive lobe of the dx2_y2 order parameter to a nega- 
tive lobe, combined with Andreev reflection by the sud- 
den TT phase shift (sign change), leads to a zero-energy 
bound statciiSiiii^ In the case of a specular [110] sur- 
face all trajectories for quasiparticles are associated with 
a sign change of the order parameter, and thus the spec- 
tral weight of the Andreev states is very large. 

In this paper we investigate theoretically the charge 
current through voltage-biased normal metal-insulating- 
barrier~d-wave superconductor (NIS) junctions subject 
to an external magnetic field directed along the c-axis. 
We extend the theory for the current fluctuations of 
conventional NIS junctions by Khlusi^ to voltage-biased 



NIS junctions with unconventional pairing correlations, 
as well as the effects of field-induced surface currents on 
the current fluctuations. We present calculations of both 
the mean charge current and the charge current fluctua- 
tions, which at zero temperature reduce to the shot noise. 
We show that the resonant nature of Andreev reflection 
via the surface bound states in a clean dx2_y2 super- 
conductor can be used to extract additional information 
from the shot noise that cannot be obtained from conduc- 
tance measurements alone. We discuss how the surface 
Andreev bound states (ABS) are broadened by impu- 
rity scattering. The impurity effect is reflected both in 
the conductance, which is related to the local density of 
states at the interface, and the shot noise, which reduces 
to 5 = 2el when impurity scattering dominates the in- 
trinsic width of the surface ABS. 

There are theoretical reasons to expect an additional, 
subdominant pairing state, e.g. with s- or dxy sym- 
metry, to be present in equilibrium under favorable 
circumstancesiSiiSiiLiSiiS The formation of surface states 
at the Fermi level, in combination with an attractive, sub- 
dominant pairing interaction favors a mixed-symmetry 
order parameter, e.g. a surface phase with dx^-y'^ ± is 
or dx^-y^ ± idxy symmetry. The sub-dominant compo- 
nent is predicted to have a phase of ±tt/2 relative to the 
dominant dx2_y2 component, and thus to exhibit spon- 
taneously broken time-reversal symmetry (T-symmetry) . 
The internal phase of the order parameter leads to a shift 
of the bound state energies below the Fermi level, thus 
lowering the surface free energy and generating a spon- 
taneous surface current ^iii 

The splitting of the Andreev states also leads to the 
prediction that the ZBCP should spontaneously split as 
a function of voltagCpS Such a splitting has been observed 
in YBCO near optimal dopingASS However, in contrast 
to the general consensus about the dx2_y2 pairing sym- 
metry for the bulk phase of the cuprates, the nature of 
the surface phase, including the possibility of broken T- 
symmetry, is unsettled liSiSl In the following we also dis- 
cuss the "fingerprints" of sub-dominant pairing that may 
be observable in the shot noise. 

The theory of shot noise in mesoscopic, metallic sys- 
tems has been used in several recent experiments to gain 
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information about tunnelling in mesoscopic systems, see 
e.g. the recent review in Ref. Earlier predictions 
for shot noise in d-wave NIS junctions were published by 
Zhu and TingiS It was shown in a non-self-consistent 
calculation that for a clean d-wave superconductor the 
Andreev resonance associated with the surface Andreev 
bound states sup pre ssed the shot noise to zero. It was 
confirmed in Ref. that this effect also holds when self- 
consistency of the order parameter is taken into account. 
Recent theoretical work on noise in d-wave SIS junctions 
addresses the effects of surface disorder, interface states 
and magnetic fields on multiple Andreev reflection in this 
system, which is an important mechanism for noise in 
voltage-biased Josephson j unctions iSSiSk 

Here we report a detailed study of the shot noise, and 
specifically address physical conditions that are likely 
present at an NIS interface. The shot noise is shown 
to be particularly sensitive to the spectrum of surface 
states and to disorder. We present results for the mag- 
netic field dependence of the shot noise, and demonstrate 
the sensitivity of the results obtained for clean d-wave 
superconductors to changes in the low-energy electronic 
spectrum by equilibrium surface currents and impurity 
scattering. We focus on the [110] orientation of an NIS 
interface to the cuprate superconductor, since for this 
orientation the influence of the Andreev bound states is 
most pronounced. 

In Section ^1 we describe our model of the normal 
metal- unconventional superconductor contact, beginning 
with a brief review of the quasiclassical Green's func- 
tion technique that we use to compute observables. We 
discuss the boundary conditions for the non-equilibrium 
propagators, coherence amplitudes and distribution func- 
tions, and express these boundary conditions in terms 
of generalized scattering amplitudes. Explicit solutions 
are used to obtain results for the transport current and 
spectral density for current noise under nonequilibrium 
steady-state conditions. In Section ITTTl we present the re- 
sults for the shot noise in voltage-biased NIS junctions 
with d-wave pairing symmetry for junctions with disor- 
der near the interface. We discuss the effects of magnetic 
fields and screening currents and surface phase transi- 
tions on the noise spectrum. Throughout the paper we 
use units for which h = ks = 1, and we choose the sign 
convention e = — |e|. 



II. THEORY AND INTERFACE MODEL 

To compute the average current and fluctuations of 
the current we use the quasiclassical Green's function 
method)2L2& and the Keldysh formalism to calculate non- 
equilibrium properties. The relevant information about 
the spectrum of current-carrying states and their distri- 
bution functions are contained in a set of non-equilibrium 
matrix Green's function: the retarded (R) , advanced (A) 
and Keldysh (K) propagators, g^'"^'^(py, R; e, t), which 
are 4x4 matrix propagators in the combined spin and 



particle-hole space (Nambu space), that depend on the 
Fermi momentum, p^, the excitation energy, e, and space 
and time coordinates, R and t. The quasiclassical prop- 
agators are related to the full Nambu propagators: 

G^(x,x') = -t&{t-t'){{'f{x),n^')}), 

G^(x,x') = -z([vl/(x), (1) 

where the Nambu field operators, ^'(x) = 

(V'T(a;),7/'i(a;),V|(x),Vx(a;))*', and = ^^x) 

incorporate particle-hole coherence of the superconduct- 
ing state. We use the notation defined in Ref. for 
the two-point functions where 9(<) is the Heavyside 
function, {A, B} = AB + BA, and [A, B] = AB - BA. 
We also use the short-hand notation x = (x, t) . 

A compact formulation of the non-equilibrium equa- 
tions is obtained in the Keldysh formulation in which 
the set Nambu-matrix propagators, Q^'^^^ ^ are grouped 
into a 2 X 2 Keldysh matrix, 

G(p,R;e,t) = J dre-^(P— ) ) . (2) 

It is most convenient to work in terms of the center-of- 
mass and relative coordinates, R = {x + x')/2 = (R, 
and r = x — x' = (r, t), and Fourier transform with 
the relative space and time coordinates. The quasiclas- 
sical propagators are then defined in terms of an inte- 
gration of the full Keldysh-Nambu matrix propagator, 
G{p, R; e, t), over an energy shell that is small compared 
with the Fermi energy, \vf{p — pf)\ < Sc <C -E/, 

g(p/,e;R,t) = i / d^p f3G(p, e; R, i) . (3) 

The quasiclassical propagator is defined by dividing out 
the weight of the quasiparticle pole in the spectral func- 
tion, a, and by convention prc-multiplying by the matrix, 
73 = fal. We denote a Nambu matrix with a 'hat', while 
Keldysh matrices are denoted with a 'check'. Thus, fs 
is the third Pauli matrix in the particle-hole sector of 
Nambu space, and 1 is the identity Keldysh matrix. 

For pure spin-singlet pairing considered here the qua- 
siclassical propagators, g^^^^^ ^ may be parameterized in 
particle-hole space by scalar amplitudes for the diagonal 
(quasiparticle) and off-diagonal (Cooper pair) propaga- 
tors, 

where X = (i?,A, if). We consider the case where the 
diamagnetic coupling of the charge currents to the mag- 
netic field dominates the paramagnetic Zeeman coupling, 
in which case the spin degrees of freedom are inert. The 
components of the quasiclassical propagators are then de- 
fined in terms of spin scalar diagonal propagators, g-^ 
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and , and spin-singlct off-diagonal propagators, 

and . These components are not all independent, but 
are related by symmetries that follow from the fermion 
anti-commutation relations^ 

For calculating the low- frequency (w ^ A, eV) con- 
ductance and noise in NIS tunnel junctions we need 
only time-independent propagators. The steady-state 
nonequilibrium quasiclassical Keldysh-matrix propagator 
obeys a matrix transport equation. 



S, 5] 



V,9 = 0, 



(5) 



where v = e$i is the electrostatic potential and S = 
A-|-Si„p represents the order parameter and impurity self 
energy. The transport equation is supplemented by the 
normalization condition, = — 7r^i)22iffi and by bound- 
ary conditions connecting the propagators at the inter- 
face. When there is no reason for confusion, we do not 
display the dependence of g and S on (p/,R;e). Equa- 
tion © represents coupled equations for the retarded, 
advanced, 



+ ivf ■ V5 



0, 



(6) 



and Keldysh propagators. 



where H 



R,A 



zv/.Vr=0, (7) 



eT3 — f — S"'"" is defined in terms of the 
excitation energy, e, the coupling to external fields, v, 

" R A . . 

and the self-energies, S ' . Similarly, the normalization 
condition expands to 



g^-^g"'^^~n^i 



g^g^-g^g^ = Q. 



(8) 



The retarded and advanced functions determine the 
particle-hole coherence functions and spectral properties 
at the NIS interface, while the Keldysh function contains 
additional information on the non-equilibrium distribu- 
tion of these states. 

The computation of g involves solving the quasiclassi- 
cal transport equations for the normal metal-insulating 
barrier-superconductor system together with a set of self- 
consistency equations for the impurity and pairing self- 
energies and boundary conditions in the bulk reservoirs 
and at the interface. 



A. Pairing Model 

The pairing correlations arc described by the pairing 
self-energy, A. In the leading order (weak-coupling) ap- 
proximation, the Keldysh component vanishes and the 
retarded and advanced self-energies are independent of 
energy. The resulting self-consistency condition is the 
BCS gap equation. 



dt 



A(p;,R) = (A(p/,p^.) / — r(p'pR;e))^, . (9) 



where /^(p/,R;e) is the off-diagonal quasiclassical 
Keldysh propagator. Wc consider spin-singlct supercon- 
ductivity derived from a pairing interaction of the form. 



(10) 



where the sum is over the 'relevant' irreducible repre- 
sentations of the crystal point group, D^h] oc € {A^^{s — 
wave)i3ig((ij.2_,y2— wave), B2g^{dxy—v/ave), ^2g(s— wave)}, 
and Aq, and ?/q(p/) are the corresponding eigenvalues and 
eigenfunctions for pairing in channel a. The 'relevant' 
channels are defined by the dominant attractive eigen- 
values for each irreducible representation obtained from 
solutions of the linearized gap equation with the micro- 
scopic pairing interaction (c.f. Ref. HqI) . The mechanism 
for pairing in the cuprates is not a solved problem, and 
indeed there may be more than one mechanism at work in 
the cuprates account for the wide range of superconduct- 
ing properties as a function of doping and disorder. For 
example, a relatively simple two-channel model based on 
electronic coupling to anti-fcrromagnctically correlated 
spin-excitations and to phonons leads to dominant d^2_y2 
pairing over a wide range of doping, but with significant 
sub-dominant pairing interactions in all other symmetry 
channels^ These sub-dominant pairing channels are pre- 
dicted to play an important role in the local electronic 
structure of surface superconducting state near an insu- 
lating barrier or other intcrface^^i^ 

In this paper we consider the signatures of sub- 
dominant pairing in the shot noise. For this purpose we 
assume the dominant pairing channel has d^2_y2 sym- 
metry, and consider sub-dominant pairing in the s- or 
do,y pairing channels, i.e. Abi, > Assg A^ig- The 
channel may also have an attractive eigenvalue for spin- 
fluctuation dominant pairing, but this order parameter 
is suppressed on both [110] and [100] boundaries, and is 
particularly sensitive to surface disorder, so we do not 
consider this sub-dominant channel for NIS junctions. 

Below the superconducting transition temperature, Tc, 
the order parameter amplitude is proportional to the Big 
basis function, rjBig = V^iPx ~ Py)- However, even a 
small attractive sub-dominant eigenvalue, Abj^ or A^^^ 
can, at low temperature generate a transition to a state 
with a mixed symmetry, with an order parameter that ac- 
quires an additional component proportional to the cor- 
responding eigenfunction, rys^^ — \/2pxPy or rjAig = Im^ 
Thus, in general we write the order parameter as 



A(p/,R) = ^A„(R)r,„(p/). 



(11) 



The gap equation separates into three self-consistency 
equations for each relevant pairing channel, 

A4R) = Xo.{vc.{Pf) £^ |^/^(p/,R;e))^^ , (12) 

The solution of the transport equation and boundary 
conditions lead to coupling between the components. 
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FIG. 1: Mixed-symmetry order parameters at a [110] specu- 
lar surface. The relative phase between the d^2_y2 and sub- 
dominant dxy or s components is ±7r/2. The insets show 
the surface bound state spectrum for these mixed-symmetry 
phases. The calculations were carried out in the clean limit 
with F = O.OOlTc and a — 10~* (see text) at a temperature 

T = omn. 



Aq,(R). The cutoff Ec and pairing interaction, Aq, are 
eliminated in favor of the instability temperatures, Tea, 
using the solution of the linearized gap equation, = 
ln(T/Tca) + /(rfe/2e)tanh(e/2r). The overall phase of 
the order parameter (|ll|l can be eliminated for an NIS 
system, but the relative phases between the different 
components that remain are determined by the minimum 
of the free energy. At sufficiently low temperature, the 
lowest energy state near a [110] surface is always a mixed 
symmetry phase with spontaneously broken T symme- 
try, in which the sub-dominant order parameter acquires 
a finite value with a relative phase of ±7r/2i2iii Conse- 
quently, we consider three possible order parameters: 1) 
pure dj.2_y2^ 2) d^^-y^ + is, and 3) dx2_y2 -\- idxy Cases 
2 and 3 are illustrated in Fig. where the pairing 

interaction parameters are chosen so that the order pa- 
rameter in the bulk region is always pure dx2_y2, and 
the subdominant components are stable near the surface 
within a layer of the order of a few coherence lengths. 



B. Magnetic Field and Screening Currents 

One of the key features of the ABS interpretation of the 
zero-bias conductance peaks observed in ab-plane tun- 
nelling spectroscopy is the splitting of the ZBCP for low 
magnetic fieldsA^SiS The energy of the ABS is shifted 
away from the Fermi level by screening currents. The 
origin of this effect is the coupling of the quasiparticle 
current to the superflow field generated in response to 
the magnetic field. 



where the condensate flow field is given by the gauge- 
invariant gradient of the phase, p,, = ^(Vi?— f A), where 
A is the vector potential. We include this coupling here 
in order to investigate the sensitivity of the noise spec- 
trum to the spectral shift of the surface ABS. Indeed as 
we show below the shot noise is particularly sensitive to 
the Doppler shift of the zero energy surface states. 

The condensate flow fleld, p^, is calculated by solv- 
ing Maxwell's equation, self-consistently with the surface 
current, supplemented with the boundary conditions for 
magnetic fleld, i? for the Meissner state far from the 
surface, and B H^^^ at the surface. For strong type II 
superconductors, such as the cuprates, with a magnetic 
penetration length A ^ ^o, the solution of Maxwell's 
equation to leading order in the small parameter ^o/A is 
Ps{z) = PsqC'^/^ with 



PsoVf _ H,. 



Ho 



dzjp{z) , 



(14) 



where jp{z) is the paramagnetic part of the current flow- 
ing parallel to the interface [computed via Eq. (|31|l be- 
low]. The paramagnetic current originates from the 
Doppler splitting of the ABS, which preferentially fa- 
vors the paramagnetic response from the bound states^ 
The field scale in Eq. is Hq = <i>o/(7r^o-^), where 

$0 = c/2|e| is the fiux quantum. 



C. Impurity self energy 

The anisotropic order parameter (Eq. is sensitive 
to disorder. We include the leading order (in l/pfi^^^ 
where l-,^^ is the mean free path) effects of disorder within 
the model of isotropic scattering of quasiparticles by im- 
purities (c.f. Ref. I33) . In this model the impurity self- 
energy is given by the quasiparticle-impurity i-matrix 
and the average impurity concentration, n.^p, 



S,„p(p/, R; e) = n.„pi(p/, p/, R; e) , 



(15) 



where f(p/,pj,R;e) satisfies a Bethe-Salpeter equation 

for repeated scattering of quasiparticles by impurities^ 
For isotropic impurity scattering defined by a scattering 
amplitude, uq, the t-matrix equations for the retarded 
and advanced self-energies have the solutions. 



i«'^(R;e) = 



and the Kcldysh component is given by 



(16) 



(17) 



The scattering amplitude, uq, defines the s-wave scat- 
tering phase shift, = arctan(7riVjMo)i while the im- 
purity concentration and normal-state density of states 

/{irNf). We use the 



''imp 



(13) 



define an energy scale Fq 
more common parametrization of the impurity scatter- 
ing model in terms of the dimensionless scattering cross 
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section, a = sin-t^o, and the pair breaking parame- 
ter, r = rosin^(So, or equivalently the mean-free path, 
£,„,=^;//2r. 

The impurity self-energy renormahzes the excitation 
spectrum via, e e(R) = e — E|*(e,R), where Y,^ is 

the fa-component of 2i„,p. For a pure d,j.2_y2 pairing 
state, the impurity renormahzation of the order param- 
eter vanishes by symmetry; the effects of pair-breaking 
are included through the renormalized excitation spec- 
trum. As a result the solution of the gap equation shows 
a reduction of the order parameter amplitude with in- 
creasing pair breaking parameter, F. For mixed symme- 
try pairing, e.g. at the surface, and in the presence of 
field-induced screening currents, the pairing self-energy 
is, in general, non- vanishing and must be calculated self- 
consistently with the renormahzation of the excitation 
spectrum. 

In general the self-energy also includes electron corre- 
lation effects generated by electron-electron and electron- 
phonon interactions. In what follows we consider a sim- 
plified model for the metallic electrodes in which the 
quasiparticles are governed by an effective one-electron 
Hamiltonian with a parabolic band structure. Thus, 
the only electronic correlations included here are those 
that contribute to the effective mass, m*, and give rise 
to superconductivity. In this case, v/ = pf/m*, and 
the charge current carried by a normal quasiparticle is 
ev/ = ^pf. Both the current and the noise spectrum 



are then calculated in this effective one-electron theory, 
modified to include pairing correlations in the supercon- 
ductor, effects of screening currents on the surface excita- 
tion spectrum and impurity scattering in both electrodes. 



D. Current and Current-Current Correlations 

Physical properties, such as the local excitation spec- 
trum, current response and correlation functions are ex- 
pressed in terms of the quasiclassical Green's function. 
Here we are interested in computing the charge current 
and the mean-field fluctuations of the current for the NS 
junction. The junction current is an expectation value, 
in a nonequilibrium ensemble (p), of the Heisenberg op- 
erator for the charge current. 



J(ri,ii) 



lim 



(Vi - V2) 



2:1 2m*i 



(18) 



Fluctuations of the current are related to the current- 
current correlation function, which is defined in terms of 
the operator, 

S(r, t, t) = K(r, t, i + r) + K(r, t + r, t) , (19) 



K(ri,<i,ti) - (-^)\lim lini (Vi - V2) (Vi - V2) 

\2m*l/ ri^ri^2^a;i\ 



(20) 



il}]{x2)tp]{xi) - ■ip\{xi)il)i{x2) 'il)\{x2)tp]{xi) - ii)\{xi)il)i{x2) - j (r i , ti )j (ri , i i ) , 



where j — Tr[/3 J] is the expectation value of current op- 
erator. The observable noise is found by evaluating the 
average of this operator over the statistical ensemble. 

In what follows we consider the effects of scattering by 
point impurities (s-wave) in the superconducting elec- 
trode with d-wave pairing. For the steady-state con- 
ductance and current noise the statistical average for 
the noise reduces to the product of two-point correla- 



tion functions. Vertex corrections to the current-current 
correlator vanish in the above approximations, or orig- 
inate from quantum interference effects or coupling to 
collective modes and thus are higher order in l/pfi^^p 
or and neglected here. Thus, after integration 

over the cross section of the junction, A, we obtain the 
current noise, S{zi,ti,ti), 



S=^ (^)'rlim^ / rfViii / d^fiii (a,, -a,J(aj, -aijTr{G<(fi,r2;fi,ii)G>(ri,f2;ii,ii)}, (21) 

f2— >f 1 

f i-»ri 

I 



where G^ = ± (G^ — G^) are the particle (<) and We separate out the momentum component parallel 
hole (>) correlation functions. to the junction using the inverse Fourier transformation 
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with respect to the difference coordinate ry 



G(ri,r2) 



i2ny 



«G(zi,z2,P||,R|j). 



(22) 



and assume that the dependence on R|| is slow on the 
scale of the coherence length, i.e. locally planar; thus, we 
omit R|| . Near the junction, incident and scattered waves 
interfere on a scale given by the inverse Fermi momentum 
pj^. We then make the following ansatz^^i^^i^ which 
factors the propagator into rapidly oscillating incident 
and reflected waves with wavenumbers, ip/z, and slowly 
varying two-point envelope functions, 

f3G{zi,Z2) = —J2 ^^(^1, ^2)e^^^''(''^^-^^=) , (23) 

where z is the coordinate normal to the interface. The 
sum is over direction indices u and fj, which are +1 or 
— 1, depending on the sign of the projection of the Fermi 
momentum on the z-axis for incident or reflected waves 
(see Fig. When v appears as an index of a function 
we use a shorthand notation ± for ±1. The diagonal 
(in direction index space) functions Cuu are related to 
Shelankov's two-point quasiclassical Green's functionsi^ 
In the limits, Z2 zi±0, these components are related 
to the projection operators. 



T iC++ =P^ = l[iT 
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-in 



(24) 



The functions Ci,(_y) are the pre-factors of the rapidly 
oscillating carrier waves, ^ g^^vtzz^ These amplitudes 
are drones - slaved to the quasiclassical propagators and 
ultimately eliminated from the boundary conditions and 
observables such as the current noise. In particular, Za- 
itsev's boundary conditions^ are derived (for details see 
Ref. l35|) by eliminating the drone amplitudes from a set of 
linear relations connecting the quasiclassical projectors, 
Ci/j/ and the quasiclassical drones, C'i,(_,y). This proce- 
dure transforms the linear boundary conditions expressed 
in terms of the set, {C,y^}, into a set of non-linear bound- 
ary conditions for the quasiclassical projectors, Ci/i/, and 
consequently for the quasiclassical Green's functions. 

For the average current, the expectation value of the 
operator in Eq. lfTH|) can be expressed as 



/(Zl,il) 

Tr{G^(a;i,X2) 



lim 



8m* X2~*xi _ 

G^(xi,x2)-G^(xi,x2)]} ,(25) 



T3 



for the current flowing through the junction along the 
z-axis. When we insert the quasiclassical envelope ex- 
pansions, Eqs. (1221)- (123, the derivat ives produce a factor 
ipfz{v + ^) ~ 2ipfz6^fi, and the current takes the form 



I{zi,t) =A— lim 

Z Z2— »Zl 



P\\ 



Tr 



(27r)2 

[T3C^^{zi,Z2,Pfe,t)^ , 



de \ ^ 



(26) 



where we neglect terms where the derivative act on qua- 
siclassical Green's functions since they are down by a fac- 
tor {pf^o)~^ compared to the leading term above. Also 
note that the term G^ — drops out because spectral 
current density is odd in energy to this order in {pf^o)~^- 
When we insert the quasiclassical expressions, 
Eqs. lj221)-(I23), into the expression for the noise, Eq. lf?T)) . 
the derivatives produce a factor {ipfz)'^{i> + p,){iy + jl) ~ 
(2ipfzY5f}^5fi,v, where the indices without (with) a tilde 
belong to the first (second) Green's function. Thus, to 
quasiclassical accuracy we obtain, 

5(zi,ti,ti)=^^_lim ^i^M X (27) 

Tr{G<^(5i,zi,P||;ti,ti)T3G>^(^i,2i,P||;^i JO-Ta} , 

where A is the cross-sectional area of the junction - which 
is the result previously obtained by Khlus^^ for conven- 
tional NS jimctions. 

The spectral density of the noise is given by the Fourier 
transformjS^ S{z,t,ijj) = / dre^'^'^ S{z,t,T), which for 
general non-equilibrium conditions depends on time. 
Here we consider the d.c. limit for the voltage-biased 
junction. This steady-state limit is independent of time, 
so we drop the time argument from here on. 

The drones, as well as the quasiclassical correlation 
functions, Gi5iy(zi, zi) = ±{C^ — C^), are continuous 
and single- valued for zi = zim^ Thus, we may take the 
limit zi —fZi; the diagonal functions arc proportional to 
the corresponding quasiclassical Green's functions. 



lim C%{zi,Z2) = 



2tt 



(28) 



and we introduce the following notation for the drones. 



Z2 



^ ^ 1^^/ \ {z^ ) 

lim C^. ,(zi,Z2)^ "^""^ 



(29) 



Thus, the current is written as 



"4/1^ E -Tr(.,.f35.a, (30) 



which is the z-component of the more general quasiclas- 
sical result for the current density. 



(31) 



where Nf is the normal-state density of states at the 
Fermi level and (. . .)p^ = Ny^ J ((27p-^^^(- ■ •) denotes 
a Fermi surface average. 

Similarly, the local noise spectrum may be written as 



i/=±i 



(32) 



7 



r R R W If 1 


92. 


r R 


~R K " K 1 








Z 


gi,[Y^,?^x^x^ ] 


92- 




^ r 2 ^ X 2 , X 2 ^ 



FIG. 2: We label Green's functions by an index 1 for the 
left (normal) electrode and 2 for the right (superconducting) 
electrode, and a direction index ± which denotes the sign of 
the projection of the Fermi momentum pf on the z-axis. The 
arrows on each line indicates the direction of p/. Also shown 
is the notation for the propagators, coherence amplitudes and 
distribution functions. 

where {. . denotes a Fermi surface average restricted 
to sgn(p/ ■ z) = v. Although the noise formula in Eq. 
(|32|l depends on the drones, they do not require indepen- 
dent calculation from their equations of motion, but arc 
expressed in terms of the quasiclassical propagators (see 
below). We note that in writing down Eq. (|32|1 . the limit 
Lo ^ was taken. 



E. Interface Boundary Conditions 

In order to compute the transport current and noise 
spectrum we must solve transport equations for the qua- 
siclassical propagators, and drones, with appropriate 
boundary conditions describing the junction. Quasiclas- 
sical boundary conditions describing partially transmit- 
ting interfaces between conducting electrodes were de- 
rived by Zaitse\i^ and Kieselmann^ for non-magnetic 
junctions, and by Millis et ali^for magnetically active in- 
terfaces. Zaitsev and Kieselmann's boundary conditions 
are a set of nonlinear equations connecting the quasiclas- 
sical propagators for incoming and outgoing trajectories 
at the interface. The material parameters entering these 
boundary conditions are the reflection (TV) and transmis- 
sion {V = 1 — 7?.) probabilities for quasiparticles when 
both electrodes are in their normal-states. This formula- 
tion is valid for arbitrary transparency. 

The quasiclassical boundary conditions, when one or 
both electrodes arc superconducting, incorporate the 
effects of particlc-holc coherence of the excitations in 
the superconducting electrodes, and the presence of 
additional channels for reflection and transmission via 
branch conversion scattering between particle- and hole- 
like branches of excitations. 

A powerful method for handling the interplay between 
these coherence effects and interface scattering was pro- 
vided by Eschrig's reductionS^ of Zaitsev and Kiesel- 
mann's boundary conditions using Schclankov's projec- 
tion operators^ and the Ricatti parametrization for the 



quasiclassical propagatorsi^SiaSiai The boundary condi- 
tion is expressed in terms of coherence functions, 7^ and 
7^, and distribution functions, and . The coher- 
ence functions have a natural interpretation as local am- 
plitudes for branch conversion; 7^ for h ^ e and 7^ 
for e h. Below we express these boundary conditions 
in terms of generalized scattering amplitudes, which in 
the clean limit are directly related to well known scat- 
tering amplitudes found in scattering theory^ The con- 
siderations below are valid for general non-equilibrium 
situations, but the results presented below are limited to 
time-independent states. For non-stationary states, all 
multiplications are replaced by time-convolutions, which 
in general prevents analytic computations. 

We adopt the notation used in Ref. |3^ for the co- 
herence amplitudes and distribution functions. The la- 
belling for functions defined on incoming and outgoing 
trajectories is also indicated in Fig. |21 For the distribu- 
tion functions, the boundary conditions can be written 
as 

^Rgx^+ff^x^ + {~T^i)i^, (33) 

= ^^.-^f + i-T/^eM + T^ui^. (34) 
-r,^xf + (-T«)if + i?f,xf, (35) 

=i-T^e)^¥+T^Hi^ + RLiL (36) 
where the scattering amplitudes are deflncd as 



' ee 




'•ee 




tR 
''eh 




(37) 


^hh 


= R^i/r, 


TR 
''hh 


= D^i/d, 


TR 

'■he 


= rdA^i, 


(38) 


f« 

ee 


= R^i/r, 


''ee 


= D^i/d, 


,R 

'■eh 


= rdA^i, 


(39) 


^hh 


- Rii/r, 


,R 

i^hh 


- D^i/d. 


,R 

''he 


= rdA^i ■ 


(40) 



The right-hand sides of Eqs. H37I4U|I are defined in Eqs. 
(D1-D5) of Ref. [H. Note that all quantities depend on 
trajectory angle pf and energy e, but not on spatial coor- 
dinates R since they arc evaluated at the junction. The 
normal-state tunnel barrier transmission and reflection 
amplitudes are denoted d and r, respectively, while the 
corresponding probabilities are denoted by V ~ d^ and 

The effective transmission and reflection amplitudes, 
including Andreev scattering, are denoted by t^^ and 
r^p, while the corresponding probabilities are denoted 
as T^p and For example, rj^f^ is the amplitude for 

reflection of a hole on the left side of the junction, while 
is the amplitude for reflection of an electron on the 
right side. Similarly, t^^ is the transmission amplitude for 
an electron from the right side to the left side, including 
branch conversion into a hole. All quantities with a bar 
refer to excitations originating from the right electrode. 

The remaining amplitudes are the Andreev reflections, 
which appear via the boundary conditions for the coher- 
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ence functions, 



"he — 


ff 


= rhhliT - 


^ thhl2d, 


(41) 


^eh — 


rf 


= reelir ^ 


- teel^d. 


(42) 


fR - 
"he — 


ff 


= rl!h72 r - 


\- thhlid, 


(43) 


fR - 

"eh ^ 


rf 


= f«7f r ^ 


-teelld. 


(44) 



In the Appendix we summarize the results for the An- 
dreev reflection probabihtics and scattering probabilities 
that enter the Keldysh distribution functions in Tables 
These probabilities are expressed in terms of the 



normal-state barrier transmission and reflection proba- 
bilities, and the coherence amplitudes for particle and 
hole excitations. 

We introduce the notation. 



1 



(a| = (1, -i(Tya*) , (45) 



which is convenient for evaluating observables. For ex- 
ample, to calculate the charge current we need 



Tr{f3|a)(/3|} = 2(l + a/3*), 



(46) 



where the factor 2 comes from the spin trace. The 
Keldysh Green's functions at the junction can now be 
written in a rather compact form. 



.gf; = -27rz7Vri [xf |r£)(r£| + If fi^f )(7f |ti 
gt = -27rz7V2-i [if fi|r« )(r« |fi + |7f )(7f |] , 



-27ri7Vn 



x^\rt){rt\+X^r,\^^){l^'\ri 



5f; = -2mN^' [x^fi\f^h){r%\ri+X^\^^){^^\] , 

(47) 

where we introduced the denominators, A^; = \C,i\'^ for 
i = 1..4 with 



Ci = l + 7fTf, 

C2 = l+lMh 



^+l2rL 
x.R^R 



CA = l+l^r^h 



(48) 
(49) 



Note that all denominators of the scattering probabil- 
ities in Eqs. (jSHl)-® (see Table |llll in the Appendix) 
are cancelled by the denominators, A^;, in Eqs. H47|) : e.g. 
Rf,N,-' = (^/Ci) * (Ci/I^n = A/\Z\^ = r£. As 
a consequence, a common denominator, \Z\'^, enters all 
Keldysh propagators. 



Z = l + i?(7f72'' + 7iV) 



(50) 



This function also appears as the denominator of the re- 
tarded Green's function. Thus, the zeroes of Z deter- 
mine the local spectrum of excitations, including inter- 
face bound states, at the junction. 

We note that the scattering amplitudes defined above 
do not exactly coincide with the ones obtained in scat- 
tering theory. There are missing pre-factors, which are 



TABLE I: Scattering amplitudes at an NTS juntion. The com- 
mon denominator is Z = 1 + TZ'y^"/^ 



„R _ r{l + 7«7«) 
ee z 


„R _ T(lW'l") 
^hh - z 


= ^ 

ee z 


f:R _ r 

rhh — z 


R _22l 

he z 


R _ 2j1 

eh z 


The = 7^7« 


= 7^7^ 


,R _ d 

''ee z 


,R d 
T'hh — 'z 


fR — d 

''ee z 


tR d 
''hh — "z 


,R _ rdj" 
'^he z 


,R _ rd-y'^ 
''eh z 


tR _ rdj'^ 
''he z 


tR _ rd-y''^ 
''eh z 



hidden in the matrices |a)(Q!| in Eqs. (|47|l . and in the 
distribution functions and [e.g. in equilibrium, 
x^ = (1 — |7^P) tanh(e/2r)]. Inspection shows that our 
generalized scattering amplitudes can be interpreted as 
describing the scattering of locally defined excitations at 
the junction, while the factors coming from the matri- 
ces and distribution functions gives a spectral renormal- 
ization due to Andreev reflection along the trajectories 
leading up to (and away from) the interface. For exam- 
ple, when a charge current is computed, these renormal- 
izations can be absorbed into the scattering amplitudes, 
which then coincide with results from scattering theory. 
However, we retained the above definitions since they 
appear naturally in the boundary condition for Green's 
functions. 

The above considerations are applicable to stationary 
states of two coupled superconductors driven out of equi- 
librium. For the special case in which the left electrode, 
is in the normal state, 7f = 7f = 0. Using Eq. we 
obtain the current computed at the junction on the left 
side of the barrier 



/ = eA^;^ I de {vf, [xf (1 + - i?fj 

Thh)] )p^. 



£f(l 



Reh ~ R 



^2 C^he T^e, 



R \ 
■hh) 

^2 C^eh 



(51) 



Explicit expressions for the effective scattering ampli- 
tudes are given in TableQlfor the NIS junction. Equation 
H51|l is valid for arbitrary stationary non-equilibrium sit- 
uations, including spatially dependent coherence and dis- 
tribution functions. Current conservation is guaranteed 
for self-consistent calculations. 



F. Asymptotic Boundary Conditions 

In the reservoir regions, far from the junction, the dis- 
tribution functions take the equilibrium forms, shifted by 
the local potential. 



Fi{x -oo, e) = tanh[(e - eV)/(2T)] 
F2{x +oo,e) = tanh[e/(2r)]. 



(52) 



The hole distributions follow by symmetry F{e) = 
F(— e). We neglect processes where quasiparticles scat- 
tered at the junction are scattered back and impinge on 
the junction before they equilibrate. The above distribu- 
tion functions then serve as incoming distribution func- 
tions in the boundary condition at the junction. We shall 
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also assume that the transparency of the junction is suf- 
ficiently small, P ^ 1, that the current flowing throught 
the system (which is proportional to V) due to the ap- 
plied voltage is small. Then, to lowest order in V we can 
neglect the effect of the current on the order parameter 
and write 



/^^(p/,R;e)= [/«-/^] tanh 



e 

2T ' 



(53) 



The necessary relations are then 
1 



di 



2VTZ 

d2s = [Dgis - (1 + n)g2s] 



2vn 

Anidia = gisdis ~ g2sd2s, 



(58) 



which is the local equilibrium form for the off-diagonal 
Keldysh propagator. We note that these assumptions will 
be valid also for high-transparency point contacts and for 
wide junctions with transport primarily through a high- 
transparency pinhole, since the current in those cases are 
reduced by the small conducting area A <C tt^q, where 
Co = ""S /Tc is the superconducting coherence length. 

Under these assumptions, the interface distribution 
functions are, 



where the first two relations come from the boundary 
condition, and the last relation is derived by making use 
of the normalization condition for Green's functions [c.f. 
Eqs. (29)-(30) of Ref.|3l|. We note that we are content 
with solving for the drones on the left side. Explicit ex- 
pressions of the drones can then be written down by using 
the Green's function g written in terms of scattering am- 
plitudes in Table |l| The retarded and advanced drones 
are 



Fi 



(1 



X2 



\l''f)F2 

-il~\lY)F, 



(54) 
(55) 



where we drop the subscript 2 on coherence functions 
since they are superfluous for an NTS system. The X2 and 
X2 terms then cancel in Eq. (|51|) by the general tilde- 
symmetry, which relates any quantity q to its partner q 
as q'(p/,R;e,i) — (?(— p/, R; — e, i)*. However, the x^ 
and x^ terms cancel only at zero bias because particles 
and holes have opposite charge. 



G. Drone Green's functions and noise 

To compute the noise in Eq. (|32() we need to also com- 
pute the drone amplitudes d^(^_^y The relations connect- 
ing the drones to the quasiclassical propagators are the 
same equations used to obtain the nonlinear boundary 
condition connecting the quasiclassical propagators gi± 
and g2±^^^ Thus, we define the symmetric combina- 
tion of Green's functions on the two sides (i = 1,2) of 
the interface as 



9i. 



9i+ - 9i-, 



(56) 



"is 



-27rirfgf3, 
+2^ir«i. 



dt = 



-2^zr,tf3, 
-27rirf„*i, 



(59) 
(60) 



while Keldysh drones take the form 



B 



+ x^ r 



1 'hh' 



KtBtR* 
•^2 ^ee''/ie 



~KtR*tR 



^2 ^hh ^eh ■> 




-iOyB* 



Z.K I R. , R.*\ 
^1 Vhh + ^hh J 



(r« 



-iUyB 



+iayB* 



;^K(^R 
hh 



-if( 



R*] 
' hh J 



(61) 



where icTy is the Pauli matrix that describes spin-singlet 
pairing. 

Since the noise is expressed in terms of = g^±{g^^ 
g^) and = d^ ± {d^ — d"^), we get terms from the 
Keldysh parts which depend explicitly on the distribution 
functions, and purely spectral terms that do not contain 
any distribution functions. Thus, we separate the noise 
into two terms, S = S^~^ + , where 



and symmetric and anti-symmetric combinations of 
drones as 



dis 
dia 



+ di^+, 
di+- - di-+. 



(57) 



S 



R-A 



4[(1 

(1^ 



Rhh)] 1 



(62) 



S'<{p^;x = 0-;e) = -(xf )22 [l + R^]' - (if )'2 [l + i?f, - R^' 

- [x^f 2(T« - T«)2 - [i^f 2(T« - T« )2 

I ^K~K,\R R* , ^R„R*\2 „K K a,R tB* , „RlR*\2 , ~K ,\ R tB* „RzR*\2 

-L ;i.K^K,\ R tR* _ R tR*,2 _ z.K~Ka\ R tR* R tR*,2 , ~K aitR tR* . tRtR*,2 
+ Xi X2 ^Feh'^ee "hh^he I -^l -^2 ^Veh^eh + "hh''hh\ +^2 ^2 ^\''he''hh + '^ee^eh I ■ 
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The above results are valid for general non-equilibrium 
distribution. The distribution functions xf and if can 
always be expressed as local equilibrium distributions 
plus 'anomalous' non-equilibrium distributions. Then 
the purely spectral terms, S^~^, are cancelled exactly 
by local equilibrium terms in that do not contain a 
Fermi function. 

Equations (|51() and (|63() . combined with Tables (P) for 
the reflection and transmission probabilities are the cen- 
tral equations needed for calculating the conductance and 
noise spectrum for NIS junctions with disorder, uncon- 
ventional pairing and interface screening currents. These 
formulas are expressed in a form that is closely related 
to the wave-function-based scattering theory applicable 
to clean systems. This connection is based on the iden- 
tification between the scattering amplitudes in the wave 
function approach and the retarded Ricatti amplitude, 
7^, which, in the clean limit, reduces to the local An- 
dreev reflection amplitude, v/u. However, the Ricatti 
representation for the propagators is more general, and 
is capable of incorporating the effects of disorder and in- 
elastic scattering. In our formulation, all observables can 
then be expressed in terms of the generalized scatter- 
ing amplitudes collected in Tabled and in the tables in 
the Appendix. However, the generalized scattering am- 
plitudes are only defined in terms of the quasiclassical 
Green's functions, through the Ricatti parametrization. 
We emphasize this fact by keeping the superscript R on 
all quantities defined in terms of the retarded Green's 
function. 

In summary, to compute the conductance and noise 
spectrum in voltage-biased NIS junctions we solve 
the quasiclassical transport equations |(SJ) for g self- 
consistently with the gap equations H12|l . the t-matrix 
equations Hlt)ll7|l . and the surface coupling to the screen- 
ing currents, (Eqs. I13I14|I . We then compute the effec- 
tive reflection and transmission probabilities, and distri- 
bution functions and use Eqs. (|51|1 and H63|) to calculate 
the conductance and the noise spectrum. 



the form, 



RnSiV) = ^J^^de {vf, (6, p;)[l - i?f,(e, pf)] 

+ i?^,(e,p/)[l-i?«(e,p/)] 
+ 2i?«(e,p/)i?^,(e,P/)})p,.,; 



■z>0 



(65) 



The normal-state junction resistance is given by R^^ ~ 
2Ae^NfVf'D, where V = (cos 0I?(p/))p^..2>o is the trans- 
port barrier transparency; cos 9 = pf ■ z > is the an- 
gle of incidence measured relative to the z-axis. Note 
that in the normal-state limit for the superconducting 
electrode, the Andreev reflection probability vanishes, 
^hei^'Pf) ~^ 0' t^*^ effective e e reflection proba- 
bility reduces to i?fg(e, p/) TZ{pf), and the integrand 
of Eq. (|64() vanishes. Thus, we recover Ohm's law for 
the junction I-V characteristic. Similarly, the shot noise 
in the NIN limit is proportional to {v fzTl{p f)'D{p f))pj. , 
which for small transparency corresponds to the Schot- 
tky result S = 2el. For higher transparency the current 
noise is reduced compared to the Schottky result. 

The above expressions have in the clean limit 
the same forms as well known scattering theory 
resultsi2i22i22iSI42iMi^ Eq. ^ also agrees with cal- 
culations of spectral current densities includ ing impurity 
scattering and subdominant pairing, in Ref. IsM 

In the following we present calculations of the zero tem- 
perature conductance 81/ dV and differential shot noise, 
dS/dV, and focus on effects of magnetic fields, impurity 
scattering, and subdominant pairing. The exact angle 
dependence of the tunnelling probability is not particu- 
larly important for our purposes; so we take it to have 
the form predicted by an interface (5-function potential, 



V{pf) = Vo 



cos 



1-Vo sin^ 



(66) 



where Vq is the transparency for normal incidence. 



III. CONDUCTANCE AND DIFFERENTIAL 
SHOT NOISE 

In the following we use these results to calculate the 
conductance and noise spectrum for NIS junctions with 
d-wave superconductors. In the zero-temperature limit, 
Eq. (|51|l for the current can be written as 

eR,J{V) = el/ + 1 r de {vf, [7^(p/) 

^ J~eV (64) 
+ i?f,(6,py)-i?«(6,p;)])p^,..^^ . 

The corresponding zero-temperature shot noise, com- 
puted at z = 0^, from Eqs. (|^ and Table iQ}, takes 



A. Pure d3,2_y2-wave: effects of a magnetic field 

In zero external magnetic field, the angle resolved dif- 
ferential shot noise is suppressed to zero at zero voltage 
and has a peak at ~ I?(p/)|Ao(p/)|, where Ao(p/) = 
^Big{z oo)r]BigiPf) is the gap in the bulk. This re- 
sult is due to the resonant enhancement of Andreev re- 
flection by the surface bound state: around the bound 
state energy (e — 0) within an energy interval set by 
the bound state width wi,{pf) ~ aI?(p/)|Ao(p/)|, the 
probability of Andreev reflection is enhanced to unity 
i?^g(e = 0, p/) = 1 independently of the smallness of the 
transparency and independently of the shape of the order 
parameter near the junction. The numerical prefactor a 
is due to the reduction of the bound state width caused by 
the suppression of the order parameter near the surface. 
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eV/T 

c 



FIG. 3: Zero temperature differential shot noise as a func- 
tion of voltage in the clean limit for several different external 
magnetic fields. Inset: The field evolution of the zero- voltage 
differential shot noise. The transparency of the interface is 
T>o = 0.1. 

It was computed for small V in Ref. and can be esti- 
mated to be approximately 1/4. As the Andreev reflec- 
tion probability is enhanced to unity, the normal reflec- 
tion probabihty is reduced to zero, i?^(e = 0,p/) = 0. 
The result of zero noise at = then follows directly 
from Eq. H55|l . The suppression of dS/dV to zero at zero 
voltage for zero field is robust under angle integration 
since the zero-energy bound state is dispersionless. The 
satellite peak will be located at a voltage of the order 
{wb{pf))pf ~ 'DqTc/2tt. This noise-Jess character of the 
zero-energy bound states in a clean system was recently 
discussed in Refs. 

In an externally applied magnetic field, the screening 
currents produce a Doppler shift of the spectrum. The 
bound state resonance is shifted accordingly. The point 
of suppressed noise is then shifted to finite voltage and 
the peak in dS/dV is pushed to higher voltages linearly 
with increasing magnetic field. These characteristics of 
the field evolution of the shot noise spectrum arc shown 
in Fig. O The dispersion of the Dopplcr-shiftcd ABS 
leads to non-zero differential shot noise at all voltages. 
In particular, at zero voltage the differential shot noise 
develops with increasing magnetic field strength as shown 
in the inset of Fig. 13 



B. d^2_y2 pairing with impurity scattering: 
Andreev vs. tunnel limits 

The sensitivity of the noise to changes in the low- 
energy surface excitation spectrum implies that the re- 
sults for S{V) in clean d-wave superconductors^SiSiii^ are 
strongly modified by disorder. Here we consider the ef- 
fects of impurity scattering on the noise spectrum. 

In Fig. ^ and Fig. |31we plot the differential shot noise 



Born limit cr=10 




eV/T r/T 

c c 

FIG. 4: (a) Differential shot noise for several different pair 
breaking parameters F for scattering in the Born limit {a = 
10"*). (b) The zero- voltage value of differential shot noise and 
conductance as a function of F. The squares and diamonds 
are the numerically computed results, while the lines are a 
guide to the eye. The junction transparency is Do ~ 0.1. In 
the inset the zero-voltage value of the differential shot noise 
is plotted as a function of q (c.f. the text), in the Born limit 
[squares - same data as in the (b)] and unitary limit [circles - 
same data as in the (b) part of Fig. 



for several pair breaking parameters. F. for scattering in 
both the Born (cr ^ 1) and the unitary limits (cr ~ 1), re- 
spectively. Impurity renormalization leads to broadening 
of quasiparticle states that depends on the pair breaking 
parameter, F, and the scattering cross section, cr. 

The local self energy at the interface is different from 
that in the bulk because of the formation of surface 
bound states. In particular, the surface bound state has 
a large impurity renormalization in the Born limit, but 
is weakly modified in the unitary limit This is oppo- 
site to the state of affairs in the bulk, where scattering in 
the unitary limit is more detrimental to the d.j.2_y2 order 
parameter, and produces a low-energy impurity band in 
the density of states. 

Impurity broadening of the surface ABS reduces the 
resonant transmission in the Andreev channel, but opens 
up the single-particle tunnelling channel. In Fig. we 
plot the reficction and transmission probabilities for scat- 
tering in the unitary limit for several values of the pair 
breaking parameter. With increasing F, the reduction of 
the Andreev reflection probability R^^ (Fig.int) is accom- 
panied by an increase of the normal reflection probability 
i?^ (Fig. Et) and an increase of the transmission prob- 
abilities, both transmission without branch conversion 
Tg^(l — |7^P) (Fig. Et) and transmission with branch 
conversion 

Thei^ - h'^f) (Fig. Hi). In particular, the 
transmission probabilities acquire a resonance form, sim- 
ilar to that in the Andreev reflection channel. 

We note that probability is always conserved during 
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FIG. 5: The same as in Fig.|l|but for scattering in the unitary 
limit [a = 1). 



scattering at the interface; it can be checked that 

+ Rt + Tllil - |7^f ) + T^^il ~ = 1- (67) 

The third and fourth terms, which describe single par- 
ticle tunnelling, are identically zero in the sub-gap re- 
gion in the absence of impurity scattering, but become 
increasingly important as the impurity renormalization 
increases (see Fig.|Hll. When the Andreev resonance is re- 
duced and single-particle tunnelling becomes important, 
the differential shot noise at zero voltage becomes non- 
zero, as shown in Figs. This is in line with the phe- 
nomenological discussion in Ref. Thus, we find that 
the noise-iess character of the zero-energy surface bound 
state is quickly lost when intrinsic broadening is present. 

To quantitatively assess the importance of impurity 
scattering in tunnelling, the contribution to the width of 
the bound state from impurity broadening, which we de- 
note Wi , has to be compared with the contribution set by 
the transparency of the interface, Wb, introduced in the 
previous section. The width Wi is related to the imag- 
inary part of the impurity self energy T,^ near the sur- 
face. Unfortunately, a rigorous analytic calculation of 
Wi in which the spatial dependence of the impurity self 
energy and the order parameter are taken into account 
has so far not been carried out (see, however, the scal- 
ing analysis in Ref. IZW) . We estimate the width to be 
Wi = c|Im{s|^(e ~ 0, z ~ 0^)}|, where c is numerical fac- 
tor which corrects for the spatial dependence of the self 
energies. In the limit Wi 3> wt, which we call the 'tun- 
nel limit', intrinsic broadening is large and only single- 
particle tunnelling is important. Andreev reflection can 
then be neglected and the shot noise for low transparency 
reduces to the Schottky form S = 2e/, and does not con- 
tain any new information that can not be extracted from 
the current. On the other hand, in the limit m; ^ Wb, 
which we call the 'Andreev limit', impurity broadening 
is negligible, single particle tunnelling is suppressed and 
Andreev reflection is resonant. In this limit the shot noise 




FIG. 6: Scattering probabilities for different pair breaking 
strengths for tunnelling ranging from the the Andreev limit 
(black curves) to the tunnel limit (brown curves). The param- 
eters corresponds to unitary scattering in Fig. |SJa), and all 
probabilities were computed at an incidence angle of 45° rel- 
ative to the interface normal. The resonance in the Andreev 
reflection probability 7?he(e) at zero energy (due to the bound 
state) is broadened by impurity scattered and suppressed in 
the tunnel limit. The resonance width in the Andreev limit is 
set by the transparency of the interface. Note that the sum 
of all probabilities is always equal to one, c.f. Eq. 1671 1. 



is non-trivial. In Figs. I4I5I the crossover between these 
two regimes is displayed for impurity scattering in the 
Born and unitary limits. It is clear that the impurity 
renormalization near the surface is much larger in the 
Born limit compared to the unitary limit: the crossover 
appears for T/Tc ^ 10"'^ in the Born limit, which is two 
orders of magnitudes smaller than in the unitary limit. 
However, if we plot RndS / d{eV) at F = as a function 
oi q = Wi/{wb{pf))pf, with the numerically computed 
Im{sf-(e ~ 0,z = 0+) and an estimate c = 1/3 in the 
Born limit and c = 3 in the unitary limit, we find that 
the crossover appears near q ~ 1 in both limits, see inset 
of Fig.Hb). 



d^2_y2 + is and d^2_ 



2 -I- idxy symmetries 



Finally, we consider the signatures of a surface phase 
transition from an inhomogeneous d^2_y2 surface phase, 
to a surface state with mixed symmetry: d^^^y^ + is or 
dx2-y2 + idxy In the clean limit the noise spectrum is 
sensitive to the change of the surface excitation spectrum 
induced by the sub-dominant pairing channel. When a 
complex order parameter develops near the surface, time- 
reversed partners of the two-fold degenerate zero-energy 
bound states are shifted in opposite directions from the 
Fermi level. The positive energy bound state spectra 
for these mixed-symmetry phases are shown in Fig. (|T|). 
A surface current, and an associated spontaneous mag- 
netic field, are generated. This symmetry breaking can 
be detected in the conductance as a spontaneous splitting 
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FIG. 7: (a)-(b) Conductance and (c)-(d) difTerential shot noise 
for an order parameter with an s subdominant component 
(left column) and a d^y subdominant component (right col- 
umn) for several different interaction strengths. The barrier 
transparency is X>o ~ 0.1, and the system is in the Andreev 
limit (r = O.OOlTc, a = 10"*). 



FIG. 8: (a) Magnetic field dependence of differential shot 
noise for a d^2_y2 + is order parameter with Tc2 = O.OSTc. 
Here is F = O.OOlTc and a = 10"*. (b) Differential shot noise 
as a function of pair breaking parameter F in the Born limit 
[a = 10--*). Here is Tc2 = 0.2T^. In both cases is Vq = 0.1. 



of the zero-bias conductance peak, or as a spontaneous 
magnetic signal from the surface. 

In the shot noise, for a clean system, we thus expect 
the point of vanishing noise to disperse with angle of in- 
cidence, in a similar way the Dopplcr shift changes the 
shot noise spectrum in an applied field. In addition to 
the dispersion of the bound states with angle, there is an 
additional mechanism of dispersion due to electron-hole 
de-phasing, which appears at non-zero energies when the 
order parameter has a spatial dependence (pairbreaking 
suppression near the surface). These mechanisms of dis- 
persion of the bound state conspire to wash out much of 
the structure one might otherwise expect to observe in 
the shot noise near zero voltage. Nevertheless, for a clean 
superconductor, characteristic differences can be seen in 
dS/dV for the cases of sub-dominant s- and dxy order 
parameters. 

In Fig. Owe plot the conductance and differential shot 
noise for several different interaction strengths, in both 
the s and dxy sub-dominant channels. There is not differ- 
ence in the signatures of the two different subdominant 
components in the conductance: in both the s and d^y 
cases the zero-bias conductance peak is split and appears 
at a finite voltage related to the size of the subdominant 
order parameter. On the other hand, in the shot noise, 
there is a double peak structure in both cases, with the 
high-voltage peak bigger than the low-voltage peak for 
the s-wave case, but with a reversal in spectral weight 
between low- and high-voltage peaks for the dxy case. 
This reversal reflects the difference in the dispersion of 
the bound states for the two different pairing channels, 
which affects the point of suppressed shot noise (as well 
as the associated peak in dS/dV). The shift to finite 
voltage is larger and disperses less in the s-wave case, 
compared to the dxy case. 



Application of a magnetic field introduces additional 
dispersion, and the angle integration leads to a reduction 
of the structures. For both types of order parameters, 
the double peak evolves with increasing magnetic field 
strength into a single peak on a scale H/Hq set by the 
size of the subdominant gap, see Fig.|Sfa). 

As in the pure d-wave case discussed in the previous 
section, impurity scattering broaden the bound state res- 
onance and reduces the structure in the shot noise. Thus, 
with increasing pair breaking parameter the two peaks in 
the shot noise merge into a single peak and, in the tunnel 
limit, the differential shot noise reduces to the conduc- 
tance, see Fig.lHJb). 



IV. DISCUSSION AND SUMMARY 

In conclusion, shot noise can be a useful tool to ex- 
tract detailed information about properties of junctions 
between normal metals and unconventional superconduc- 
tors. However, a necessary condition is that the system 
is in the Andreev limit, Wi ^ Wb, as shown in Figs. 01- 
IHl This is a rather restrictive condition at present, since 
in most experiments the zero-bias conductance peak is 
broadened by disorder, sec however Rcf. 

In several recent experiments^i^ the density of states 
around single impurities or inhomogeneities were mapped 
out by c-axis STM spectroscopy. The results are dis- 
cussed in terms of low energy states bound to a single 
impurity scattering in the unitary limit, in line with the- 
oretical works in Refs. lilsnllsiL Thus, if impurities in the 
high- Tc superconductors are indeed scattering in the uni- 
tary limit, we expect that the surface bound states will 
not be particularly broadened, and tunnelling in the An- 
dreev limit should be possible to achieve experimentally 
in a clean sample. The mean free path corresponding to 
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g ^ 1 in the unitary limit (for the parameters in Fig. |SJ) 
is estimated to be of the order of tenths of coherence 
lengths, which is achievable experimentally. 

There are other sources of broadening of the sur- 
face/interface bound states that were not considered in 
this paper. In particular, it is clear that surface rough- 
ness will drive the system towards the tunnel limit, be- 
cause non-specular scattering of quasiparticles to the 
nodes of the order parameter broadens the bound states 
just as impurity scattering does. Therefore, to extract 
information from shot noise it will be important to have 
a specularly reflecting junction, or tunnel from an STM 
tip directly into the a6-plane of a specular portion of a 
superconductor surface. 

Under these circumstances, information about the su- 
perconductor properties can be deduced via the partic- 
ular properties of the zero energy surface bound states. 
The shot noise in a purely (i^2_j,2-wave superconductor 
is suppressed around low voltage and approaches zero in 
the clean limit. The characteristic magnetic field depen- 
dence shown in Fig.Ocan then be used to test the theory. 
The zero- voltage shot noise level changes according to the 
inset of Fig. 13 and the satellite peak is linearly pushed 
out with increasing magnetic field strength. The double 
peak structure shown in Fig. |7| serves as a fingerprint of 



the symmetry of the subdominant pairing channel. 
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Appendix 

In this appendix we tabulate the Andreev reflection 
probabilities expressed in terms of the barrier reflection 
and transmission probabilities and coherence amplitudes 
for incoming trajectories in Table UTI and in Table ITTTI we 
summarize the transmission and reflection probabilities 
that enter the boundary conditions for the distribution 
functions in Eqs. H33I36|I . 
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